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Motivated by recent astrophysical observations of quasar absorption systems, we formulate a 
simple theory where the electron to proton mass ratio p = m e /m p is allowed to vary in space-time. 

In such a minimal theory only the electron mass varies, with a and m p kept constant. We find 
that changes in p will be driven by the electronic energy density after the electron mass threshold is 
crossed. Particle production in this scenario is negligible. The cosmological constraints imposed by 
recent astronomical observations are very weak, due to the low mass density in electrons. Unlike in 
similar theories for spacetime variation of the fine structure constant, the observational constraints 
on variations in p imposed by the weak equivalence principle are much more stringent constraints 
than those from quasar spectra. Any time-variation in the electron-proton mass ratio must be less 
than one part in 10 9 since redshifts z « l.This is more than one thousand times smaller than current 
spectroscopic sensitivities can achieve. Astronomically observable variations in the electron-proton 
must therefore arise directly from effects induced by varying fine structure ’constant’ or by processes 
associated with internal proton structure. We also place a new upper bound of 2 x 10~ 8 on any 
large-scale spatial variation of p that is compatible with the isotropy of the microwave background 
radiation. 


I. INTRODUCTION 

This work is motivated by recent observational con¬ 
straints 0, a i on variations of the electron-proton 
mass ratio p = m e /m p which take advantage of new 
high-precision spectra to complement previous investiga¬ 
tions of p in combination with the fine structure con¬ 
stant, a, and the proton g-factor, g p , 111 . It 
is important to relate these constraints to independent 
investigations on varying a 011 The data re¬ 
ported in ref 1 from 21 cm and quasar UV spectra 
are drawn from absorption spectra in the redshift range 
0.24 < 2 < 2.04 and lead to an observational constraint 
on shifts in x = g p a 2 p of Ax/x = 0.35 ± 1.09 x 10 -5 . 
As discussed in detail in ref JJ, using data on varia¬ 
tions permitted in a by Murphy et al 7j and Chand 
et al J8j this results in constraints on variations in p of 
A p/p = 1.5 ± 1.1 x lO" 5 and A p/p = 0.5 ± 1.1 x 10“ 5 , 
respectively, since we expect that the variations in the 
g-factor A gjg p « -0.1A(m q / AQ C D)/m q /A q C d will be 
negligible [l0|. The question we ask in this paper is: how 
should theorists react should a varying p be observed but 
not a varying alpha? Could a simple framework accom¬ 
modate such a situation? 

Theoretically there is a wide range of possible connec¬ 
tions between a and p. In previously considered models 
variations in a lead to variations in the electron mass (via 
the electron self-energy) and in the proton mass (via the 
electrostatic energy contained inside a proton). These are 
model dependent and in general quite complex to work 
out (Ill, (lj, M • However it certainly seems the case that 
a varying alpha entails a varying p. 

However, superficially it should also possible to have a 
variation in p without a variation in a. The idea is simply 
to induce such variations directly in the mass parameter 
for the electron. At least at tree level, such variations 
would have no effect upon the mass of the proton, or 


the value of a. This constitutes the simplest varying 
p theory and in this paper we investigate some of its 
cosmological consequences together with the constraints 
that can be placed by considering the allowed level of 
violations of the weak equivalence theory. The analogous 
theory for variation of a developed in ref. 0 produced 
observable effects on quasar spectra at high redshift {z ~ 
1) with accompanying violations of the weak equivalence 
principle of order 10~ 13 , an order of magnitude below 
observational bounds. In contrast, if a varying p exists at 
the astronomically observable level it result in violations 
of the weak equivalence principle that are unacceptably 
large by a factor of more than 10 3 . 


II. THE FORMALISM 

We shall use metric convention and set 

H = c = 1 throughout. Consider the standard Dirac 
lagrangian 

Cy = — to’F’F (1) 

and let the electron mass be controlled by a “dilaton” 
field (/> defined by m = mo exp (f>, where mo is the cur¬ 
rent electron mass (so that <f> = 0 today). The minimal 
dynamics for <fi may be set by the kinetic lagrangian 

( 2 ) 

where w is a coupling constant. This is the minimal the¬ 
ory; we can add a mass, a potential, or any other com¬ 
plication if required. 

From this lagrangian we obtain a Dirac equation with 
variable electron (and positron) mass 

(i — m)A/ = 0. 


( 3 ) 
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The dynamical equation for the logarithm of the mass 
(</> = ln(m/mo)) is the driven wave equation 

d 2 p=- — 4/4/ (4) 

iO 

Performing a simple calculation we can derive the follow¬ 
ing result concerning the driving term: the macroscopi- 
cally averaged value of 4/4/ is negligible in the relativistic 
regime and is given approximately by the electron and 
positron number density ( n e + n p ) in the non-relativistic 
regime. Thus, unlike the situation in the simplest theo¬ 
ries for variations in cc, Ena, variations in /i do occur 
in the radiation era and start as soon as the universe 
cools down below the electron rest-mass threshold. 

The calculation is straightforward E3- We want to 
compute the macroscopically averaged value of 4/4/: 

<**)> = YKtJ d 3 xdtw . ( 5 ) 

From the standard expansion 

*(®)= / #p^(a s (pK(p)e- i ^+5t(pK(p)e i ^) 

J S =1 

( 6 ) 

we arrive at 

J d 3 x\M = j d 3 p^^( 4 (p)a r (p) - 6 r (p)^(p)) 

rs 

+ ^2(a r (p)b s (-p)v r {p)u s (-p))e~ 2tEt .( 7) 

rs 

The last two terms average to zero when integrated in 
time. The remaining terms, upon normal ordering, lead 
to 

<*4/> = ±J d 3 P^(N e (p) + N p (p)) (8) 

Thus, in the relativistic regime E m and the driving 
term is negligible. In the non-relativistic regime, on the 
other hand, m reduces to n e +n p . 

This theory conserves lepton number. The term e^4/4Z 
leads to Feynman diagrams involving <p and the 4/, but if 
electrons are created then an equal number of positrons 
must also be created. In general, a varying mass leads 
to particle production (see [n|; and compare to the case 
of pair production in an electric field 1151 1. In the non- 
relativistic regime, however, this is negligible as shown 
below. 


III. COSMOLOGICAL SOLUTIONS 

We assume that the variations in p that drive varia¬ 
tions in p are small in the sense that they do not produce 


significant contributions to the Friedmann-like equation 
governing the dynamics of the expansion scale factor of 
the universe, a(t). This assumption will be confirmed 
by the results obtained. In a cosmological setting the 
equation for <p is 

•• CL • 777 , 

<p-\-3 —<p ~ - (n e + n p )Q(m — T) (9) 

CL LU 

where 0 is the theta function and a(t) will be the scale 
factor for a spatially flat Friedmann universe containing 
radiation, matter and quintessence. The limiting case of 
no mass variations corresponds to the u —> oo limit. We 
first try to constrain this parameter cosmologically. We 
write: n e ss fnb = / 77 s, where ul, rib and s are the lepton 
number, baryon number and entropy densities, / ~ 0.95 
under standard assumptions and p = rib/s must be of 
order IO -10 for Big Bang nucleosynthesis to produce the 
observed light element abundances. 

Lepton number (and charge) conservation in this the¬ 
ory implies that til oc 1/a 3 . This imposes the require¬ 
ment that 

h e + 3— n e = P(t ) = h p + 3—n p (10) 

where P(t) is the particle production rate density. En¬ 
ergy conservation in turn requires that 

a • a 

Pcf, + 3-(p<£ +P 0 ) = -mn L (j) = —p L - 3 ~Pl (H) 

where Ptj> = P 4 , = uxfr/2 , and (in the non-relativistic 
regime) p l = m(n e + n p ), p^ ~ 0. Combining energy and 
lepton-number conservation we therefore conclude that 
P(t) = 0 and we have n e = B/a 3 and n p = A/a 3 ; A, B 
constants. Although the electron number is conserved 
the varying mass induces a corresponding variation in the 
electron’s energy density. This is supplied (or absorbed) 
by the <j> field so that the overall energy is conserved. 

We now consider a number of approximate solutions 
with cosmological implications. In general, we can write 
© as: 

(/>a 3 j = — M exp[(/>] (12) 

with 

M = a 3 n L m 0 _ p e0 ag ^ 

U) u> 

If the mass variations are small, as observations demand, 
we may set « 1 in the right hand side of <E3 and 
integrate to obtain 

<pa 3 = -Mt + A (14) 

where A is an integration constant. In the dust era, we 
take a oc f 2 / 3 , and assume that the driven solution dom¬ 
inates [ 2 (if : hence A ss 0 , and so p = — Mt\ ln(t/fo), that 
is 

<p = ^Mtpln)! + z) (15) 
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where to is the present time and z is the redshift corre¬ 
sponding to comoving proper time t. But 


n e 


Mt 2 0 = —p c t 2 0 = 


Q b m e 


1 - 


file 


1 


6 ttG’ 


(16) 


where fn e is the heliunr-4 number fraction, and so the 
shift in p between redshift 2 and the present is given by 

T a * = 0 - Sj f ) ln(1 + (17> 


The observational constraints from redshifts of order 1 
mss therefore convert into a bound on ui the one free 
parameter defining the theory and 

^ < 10" 5 -f GU > 0.2. (18) 

This is a very poor constraint, a fact to be expected be¬ 
cause the electron density is so low. It means that u> 
can be extremely ’small’ - slightly smaller even than the 
minimal scale of E 2 lanck = G _1 . Recall that the no- 
m e -variation model has u> = oo and current astronom¬ 
ical observations at high redshift barely constrain this 
maximal variation case. Similar bounds on the combina¬ 
tion a 2 p arise by introducing an additional scalar field 
to carry variations in a as in |lj] and the bounds on the 
time variation of this combination are dominated by the 
constraints on changes in p if the two scalar fields are 
coupled only by gravity. 

In the above derivation we have neglected the accelera¬ 
tion of the universe. This will weaken the bound further, 
as acceleration stabilizes p. If we solve m in a flat 
background universe containing dust and a cosmological 
constant, A, and scale factor a 3 oc sinh 2 [tv / 3A/2] then 
the driven solution becomes 

Afflfro , A/ft^o . 

</> = </>* + —-— t coth At-——— In smh At (19) 

“AO A “AO-W 


where A = -\/3A/2 and as f —> oo we have <f> —> con¬ 
stant which we neglect as small. Similar solutions have 
been found in (2J, |25j . To a reasonable approximation <j) 
evolves as G5) until the epoch z\ = -^/7/3 — 1 ~ 0.3 
when the cosmological constant begins to dominate the 
expansion and is approximately constant thereafter. So 
the bound m will be weakened by a factor of order 0.3. 
The variation of </> is turned off in a similar way when¬ 
ever the dynamics are dominated by a fluid with equa¬ 
tion of state satisfying p + 3p < 0 (such as quintessence) 
and this includes the case of negative curvature (a = t) 
also. We see that our distance from the future asymptotic 
value, 4 >(to) — </>*, is also small. 

In the radiation epoch, where a = f 1 / 2 , we find a dif¬ 
ferent solution of 1121 ) : <j>a 3 = —M(t — t m ), where t m 
ss 100s is the time when the electron mass threshold is 
crossed at T ~ m e . This integrates to give 


where the constant C should be chosen so that at equal¬ 
ity <f> = (f) eq matches the result obtained from the dust 
solution 11 7[) . As z m z eq we have 

^ = ^ eq +2Mt 2 0 ( — (21) 

\ Z e q Z J 

(there is a small correction if 2 ss z m ~ 10 9 ). We see 
that most of the variation in p in the radiation epoch 
happens near the epoch of matter-radiation equality, so 
that the redshift of the electron mass threshold, z m , is 
numerically irrelevant. 

We can now consider the total variation in p as we 
go up in redshift deep into the radiation epoch, 2 —> oo. 
This is given by 

^=0oo = Mt 2 (hnz eq + —). ( 22 ) 

M \ ^ Zeq J 

The first term contains the variation that happened in 
the matter epoch; the second adds on the variation that 
happened in the radiation epoch. The first term domi¬ 
nates, and so most of the variation in the electron mass 
in this theory occur in the matter epoch. Thus we do not 
expect detailed consideration of primordial nucleosynthe¬ 
sis at 2 ~ 10 9 — 10 10 to constrain it further. Note that 
our approximation was based on neglecting the contribu¬ 
tion of the scalar-field energy density to the Friedmann 
equation. This requires u>cj) 2 /t~ 2 to be bounded at early 
times and this holds for our solution H3). 

This leaves us with a worrying possibility. In the above 
discussion we have assumed that before t m ~ jn” 1 there 
are no variations in m e since the driving term is ab¬ 
sent. But this need not necessarily be the case. We 
could add to the driven solution any solution to the free 
equation, cj)a 3 = —A, where A is an arbitrary constant. 
In the radiation-dominated case |^J this would add to 
(j> solution cfifree = —2 A/t 1 / 2 As long as A is chosen so 
that (f> variations are suitably small at nucleosynthesis, 
this component cannot be ruled out. Its effect, how¬ 
ever, would be dramatic. It would imply that the elec¬ 
tron becomes massless (or infinitely massive) at the Big 
Bang. However, if this solution mode is included our 
approximation of neglecting the back-reaction of the (f> 
motion on the Friedmannian dynamics of the universe 
begins to fail. The kinetic term for this mode is of order 
ij(j) 2 ss wA 2 a -6 and will dominate the radiation plasma 
(p oc a~ 4 ) in the Friedmann equation as t —> 0, leading 
to scalar-dominated background expansion with a = t 1 / 3 . 
The resulting exact general solution of d is 


m = exp[0] 


2C 2 / t \ ±G 1 

Aft ^77 [1 - (t/T) ±c ] 2 


(j) = C — 2Aftn — 
ao 


1 + T 


( 20 ) 


with C constant. With C = 2n+1, n £ Z + , T = —/3 < 0, 
we have m oc t 2n —> 0 as t —> 0. 
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IV. THE WEAK EQUIVALENCE PRINCIPLE 


Stronger bounds on this theory are possible if we con¬ 
sider the observational constraints on any local violations 
of the weak equivalence principle (WEP). All varying 
constant theories imply the existence of a ’fifth force’ but 
this need not violate the equivalence principle. Varying e 
theories, like BSBM, which match the level of variations 
in a consistent with quasar data of refs [?J predict WEP 
violations within one order of magnitude of current con¬ 
straints 0 and there is a real prospect of testing this 
prediction with the first generation of space-based tests 
of the WEP. In the varying-/! theory introduced here the 
situation turns out to be quite different. 

Defining ( t = p e / /?, test particles may be represented 

by 

£(y) =- J dr m((l — Ct) + Qe 0 ) [~g^:x^x"] * ^ 

(23) 

where over-dots are derivatives with respect to the proper 
time r. This leads to equations of motion: 

^ + ^- (1-^;,,/ ^° <24) 


which in the non-relativistic limit (with (t < 1) reduce 
to 

^|! = -V i $ + C(25) 

where $ is the gravitational potential. Linearizing 0 
around a spherically symmetric body with mass M s and 
C = 0 we have 


Cs GM S 
+ku> r 


(26) 


so that we predict an anomalous acceleration equal to 


a = 



CsCt \ 
Gujtt J 


(27) 


Violations of the WEP occur because Ct is substance de¬ 
pendent. For two test bodies with £1 and £2 the Eotvos 
parameter is: 


? _ 2 |ai — 02) _ ClCi ~ C2I 

a\ + 02 Gujtt 


(28) 


It is easy to see why the varying-/! theory is so much 
more strongly constrained by tests of the WEP than the 
BSBM varying-a theory 00- The latter is a varying- 
e theory, driven by E 2 — B 2 . It is possible that the dark 
matter is dominated by E 2 or B 2 field: this enhances 
the cosmological variations permitting variations in alpha 
at z ~ 1 of order 10 -5 without conflict with the WEP, 
as explained in 0. A similar structure for the theory 
proposed in this paper would require the dark matter to 
be dominated by electrons, that is Co ~ 1. This is not 
a possibility. Still one should bear in mind this possible 
loophole, as well as the one raised in p^ ]. 

V. COSMOLOGICAL INHOMOGENEITIES 

Another astronomical consideration is that of the non- 
uniform matter distribution in the universe. In line with 
usual practice we have considered the evolution of p in 
a spatially homogeneous universe. Since the real uni¬ 
verse exhibits significant inhomogeneity in the matter 
distribution on sub-Gpc scales we expect there will be 
inhomogeneity in p also. From eqns. (121I1 - I122I) we see 
that the fact that the dominant contributions to vari¬ 
ations in p come from changes at low redsliifts in the 
dust-dominated era. The largest effects arising from the 
effects of growing density inhomogeneities will arise at 
these times. The leading-order effect in (E3> will be in¬ 
homogeneity in M which arises directly from inhomo¬ 
geneity in p e . Once protoclusters separate out from the 
expansion of the universe and virialise the time evolu¬ 
tion of p will cease inside them while p continues to de¬ 
crease logarithmically in the back grou nd universe until 
the expansion starts to accelerate [23. This will create 
small spatial variation in the value of p between clus¬ 
ters and the background universe and we can solve o 
to determine the evolution of (f> and p inside the over- 
dense region. If we model a growing spherical density 
inhomogeneitv in the dust era by a closed Friedmann 
universe 1271 128] with greater than average density we 
obtain </> = —Mt [ln(f/f 0 ) + T(x)(t/to) 2 / 3 ] where r(x) 
is a positive function of the space coordinates that fixes 
the amplitude of the growing perturbation mode in the 
linear approximation. Hence, substituting for redshift 

(j) = -Mtg ln(l + z) — _ 

v 2 0 V ’ 1 + z 


Since ( s ~ |Ci — C2I = O(l0~ 4 ) then in order to produce 
77 < 0( 10 -12 ) we need Guj > 10 3 . Similar bounds are 
obtained by considering the motion of the Earth-Moon 
system due to their compositional differences, H3 
Under this constraint on ui the largest tolerated varia¬ 
tions in p at redshifts of order 1 from eq. o are of the 
order 

(29) 


To leading order, the spatial inhomogeneity in p is 

= -Mtim (i) 1 ' 3 = ~ Mt ° rw 

p p \to ) 1 + z 

where p is the spatially uniform value of p in the back¬ 
ground (where T = 0). Therefore we expect inhomo¬ 
geneities in p to grow in time at the same rate as inhomo¬ 
geneities in the matter density which drive them. There 


\Ap\/p < 10" 9 . 
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is therefore observational motivation to search for spatial 
variations in p and in p. Since the dominant evolution 
of the electron mass occurs in the dust-dominated era 
there appears to be little effect on primordial nucleosyn¬ 
thesis. There will be effects on the Thomson scattering 
of the microwave photons (as <jt ~ m“ 1 2 ) but they are 
also expected to be unobservably small. 

It is possible to place a strong bound on the magnitude 
of any spatial variations in p by noting that inhomogene¬ 
ity in p is driven by inhomogeneity in the density p e and 
the latter will lead to metric potential perturbations on 
the scale of the inhomogeneity which will produce tem¬ 
perature anisotropies in the microwave background on 
large angular scales {2£?j ■ If we assume that the inhomo¬ 
geneity in the electron density is approximately propor¬ 
tional to that in the matter distribution, so 

$Pe j & Pm 

Pe Pm 

then the spatial inhomogeneity in the electron-proton 
mass ratio, p, is 

Sp 0.3 5p e fL\ 2 4 0.9/3 AT 2 x 10" 5 /3 2xl0 - 

p Goo p e \t) Geo T Goo 

Here, we have used the WEP lower bound on Goo. I2HI to 
improve on the cosmological bound given in 29\ on spa¬ 
tial inhomogeneity in the electron-proton mass ratio on 
large cosmological scales that have not undergone gravi¬ 
tational binding and non-linear evolution. 

VI. CONCLUSIONS 

In summary, we have formulated the first theory which 
describes the cosmological evolution of the electron- 


proton mass ratio, p. Although ’limits’ exist on the time 
variation this quantity in the literature and the conse¬ 
quences of altering its constant value are well known 
[3dL I/iTl .'Uj- 33, Hi there has been no self-consistent theo¬ 
retical description of the time-evolution of p that is con¬ 
sistent with the conservation of energy and momentum 
and cosmology. In this paper we have presented the first 
(and simplest) such theory. It can be related to dilaton 
couplings beyond tree-level found in string theory l35j , 
to the covariant varying-c model [Hj, or to the standard 
model with varying parameters (as long as the variations 
envisaged here are due to a dynamic Yukawa coupling 
rather than a shift in the Higgs VEV [UJ). The theory 
has appealing features which result in it being relatively 
insensitive to uncertainties in defining parameters and 
we have found that it is remarkably weakly constrained 
by the high-precision astronomical data provided by ob¬ 
servations of quasar s pec tra. However, unlike the case of 
theories of varying a, Il4| . the theory is significantly more 
strongly constrained by the data supporting the weak 
equivalence principle. The deviations from composition- 
independent freefall under gravity place constraints upon 
varying p which are about 1000 times stronger than the 
current quasar bounds and indicate that direct effects 
from varying p on atomic energy levels (in particular to 
the observable combination a 2 7 8 p) at redshifts 0(1) will be 
unobservably small (at most O(10~ 9 )) if the variations in 
p are driven by scalar-field variation of the electron mass 
as assumed here. Thus if observational evidence arises for 
cosmological time variations in p we expect complemen¬ 
tary variations in a to exist from which the p variations 
arise from radiative corrections or electrostatic effects in¬ 
ternal to the proton. 
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